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LONG-TERM  GOALS 

The  long  term  objective  of  our  research  for  the  “High  Resolution  Air-Sea  Interaction”  (HIRES) 
Departmental  Research  Initiative  (DRI)  is  to  identify  the  couplings  between  large  wave  events, 
winds,  and  currents  in  the  surface  layer  of  the  marine  boundary  layers.  Turbulence  resolving 
large  eddy  simulations  (LESs)  and  direct  numerical  simulations  (DNSs)  of  the  marine  atmo¬ 
spheric  boundary  layer  (MABL)  in  the  presence  of  time  and  space  varying  wave  fields  will  be 
the  main  tools  used  to  elucidate  wind-wave-current  interactions.  A  suite  of  turbulence  simula¬ 
tions  over  realistic  seas  using  idealized  and  observed  pressure  gradients  will  be  carried  out  to 
compliment  the  field  observations  collected  in  moderate  to  high  winds.  The  database  of  simu¬ 
lations  will  be  used  to  generate  statistical  moments,  interrogated  for  coherent  structures,  and 
ultimately  used  to  compare  with  HIRES  observations. 


OBJECTIVES 

Our  goals  are:  1)  participate  in  the  planning  process  for  the  HIRES  field  campaign;  and  2)  con¬ 
struct  and  exercise  an  LES  code  applicable  to  the  HIRES  high  wind  regime. 


APPROACH 

We  are  investigating  interactions  between  the  MABL  and  the  connecting  air-sea  interface  pri¬ 
marily  using  LES.  The  waves  are  externally  imposed  based  on  well  established  empirical  wave 
spectra.  Ultimately  they  will  be  provided  by  direct  observations  of  the  sea  surface  from  the 
HIRES  field  campaign.  The  main  technical  advance  is  the  development  of  a  computational  tool 
that  allows  for  nearly  arbitrary  3-D  wave  fields,  i.e.,  the  sea  surface  elevation  h  =  h(x,y,t )  as  a 
surface  boundary  condition.  The  computational  method  allows  for  time  and  space  varying  sur¬ 
face  conditions  over  a  range  of  wave  scales  (9(10)m  or  larger  as  shown  below. 


Report  Documentation  Page 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 


1.  REPORT  DATE 

30  SEP  2012 


2.  REPORT  TYPE 


3.  DATES  COVERED 

00-00-2012  to  00-00-2012 


4.  TITLE  AND  SUBTITLE 

Turbulent  Flow  and  Large  Surface  Wave  Events  in  the  Marine  Boundary 
Layers 

6.  AUTHOR(S) 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

National  Center  for  Atmospheric  Research, Boulder, CO, 80307-3000 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 


13.  SUPPLEMENTARY  NOTES 


14.  ABSTRACT 


15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF: 


a.  REPORT 

unclassified 


b.  ABSTRACT 

unclassified 


c.  THIS  PAGE 

unclassified 


17.  LIMITATION  OF 

18.  NUMBER 

ABSTRACT 

OF  PAGES 

Same  as 

11 

Report  (SAR) 

19a.  NAME  OF 
RESPONSIBLE  PERSON 


Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


WORK  COMPLETED 


Meetings:  We  attended  the  Ocean  Sciences  Meeting  in  Salt  Lake  City,  UT  and  the  20th  Sym¬ 
posium  on  Boundary  Layers  and  Turbulence  joint  with  the  18th  Conference  on  Air-Sea  Interac¬ 
tion  in  Boston,  MA  where  we  presented  results  from  research  funded  by  HIRES  and  the  com¬ 
panion  DRI  Impact  of  Typhoons  on  the  Ocean  in  the  Pacific  (ITOP).  During  the  past  year, 
we  also  collaborated  extensively  with  numerous  investigators  who  also  presented  results  from 
current  research  at  these  meetings.  The  research  spanned  a  broad  range  of  topics  in  the  atmo¬ 
spheric  and  oceanic  boundary  layers.  A  summary  of  the  presentation  titles  and  archival  manuscripts 
produced  during  the  past  year  are  given  in  the  PUBLICATIONS  section  of  this  report.  We  also 
participated  in  the  HIRES  PI  meeting  at  Scripps  Institution  of  Oceanography  held  in  April 


2012. 


Table  1:  Simulation  properties 


Run 

Ug  (m  s"1) 

A  =  Cp/U10 

u*  (m  s  x) 

Q*  (K  m  s  x) 

L 

A 

5 

4.8 

0.124 

0.0 

oo 

AC 

5 

4.8 

0.166 

0.01 

-35.0 

ACF 

5 

- 

0.167 

0.01 

-35.6 

B 

7.5 

3.4 

0.187 

0.0 

oo 

C 

10 

2.8 

0.228 

0.0 

oo 

D 

15 

1.9 

0.338 

0.0 

oo 

E 

20 

1.5 

0.452 

0.0 

oo 

EB 

20 

- 

0.471 

0.0 

oo 

EC 

20 

1.5 

0.539 

0.01 

-1200.0 

Simulations:  During  the  course  of  the  HIRES  DRI,  we  built  a  database  of  LES  solutions  that 
examine  the  influence  of  wave  age  on  MABLs  driven  by  a  combination  of  shear  and  convective 
forcing.  Often  the  wind  driven  MABL  is  either  cooled  or  heated  at  the  surface  and  thus  it  is 
of  interest  to  consider  flow  over  waves  in  the  presence  of  atmospheric  stability.  Previously,  we 
examined  the  impact  of  stratification  using  DNS  (Sullivan  &  McWilliams,  2002)  and  LES  (Sul¬ 
livan  et  al.,  2008)  over  an  idealized  monochromatic  wave.  Table  1  lists  bulk  properties  of  the 
simulations,  viz.,  the  geostrophic  wind  Ug,  wave  age  A  =  Cp/U\ 0,  friction  velocity  zt*,  surface 
heat  flux  Q *,  and  Monin-Obukhov  stability  index  L.  In  the  definition  of  wave  age,  Cp  denotes 
the  phase  speed  of  the  peak  in  the  wave  height  spectrum  (see  description  below)  and  Uw  is  the 
reference  wind  speed  at  a  height  of  10  m  above  the  water  surface.  The  height  of  the  MABL  in 
these  simulations  is  zt  ~  500  m  and  in  the  experiments  with  unstable  surface  forcing  the  Dear- 
dorff  (1972)  convective  scale  w*  ~  0.53  ms-1.  Simulations  (without  wave  effects)  are  in  a  mixed 
shear-convective  regime  that  tend  to  develop  broad  streanrwise  rolls  (e.g.,  Moeng  &  Sullivan, 
1994).  In  addition  to  the  simulations  with  a  wavy  lower  boundary  we  also  consider  a  run  ACF 
with  a  flat  lower  boundary  as  a  reference  for  run  AC  and  a  run  EB  with  a  stationary  lower 
wavy  surface  (or  bumps)  for  comparison  with  run  E.  In  all  the  LES  experiments,  the  compu¬ 
tational  domain  is  (XL,  YL,  ZL)  =  (1200, 1200,  800)  m  and  is  discretized  using  (Nx,  Ny,  Nz)  = 
(512,512, 128)  gridpoints.  A  description  of  the  code  and  first  results  are  described  in  Sullivan  et 
al.  (2010a).  We  note  that  the  algorithmic  formulation  is  general  and  also  allows  turbulence  sirnu- 


lations  over  and  around  3D  stationary  orography  as  described  in  Sullivan  et  al.  (2010b)  1 . 

A  broadband  wave  spectrum  is  considered  in  the  numerical  experiments  with  a  wavy  lower  sur¬ 
face.  The  virtual  surface  wave  field  is  built  to  match  a  wind-wave  equilibrium  Pierson  and  Moskowitz 
(1964)  spectrum  assuming  a  surface  wind  speed  of  15  m  s-1;  hence  the  phase  speed  of  the  peak 
in  the  wave  spectrum  Cp  ~  18  m  s-1.  The  individual  wave  modes  in  the  wave  spectrum  are  as¬ 
sumed  to  have  random  phases  and  a  directional  spectrum  is  picked  to  emphasize  long  crested 
waves.  Each  wave  mode,  resolved  by  the  LES  mesh,  is  advanced  forward  in  time  using  the  lin¬ 
ear  deep  water  dispersion  relation.  At  each  stage  of  the  Runge-Kutta  time  stepping  scheme  we 
construct  a  new  grid  in  physical  space  that  fits  the  wavy  surface.  At  a  particular  (x,  y )  gridpoint 
the  stretched  vertical  grid  spans  h(x,y,t )  <  z  <  Zl  where  h(x,y,t )  is  the  instantaneous  wave 
height  and  Zl  is  a  fixed  height  above  the  water  surface  (independent  of  (x,y,t)).  An  x  —  y  view 
of  the  3D  surface  wave  field  at  a  particular  instant  in  time  is  given  in  figure  1  and  an  oblique 
cutout  view  of  the  computational  grid  above  the  waves  is  depicted  in  figure  2.  Additional  details 
of  the  LES  experimental  design  are  provided  in  prior  year  progress  reports. 


RESULTS 

In  previous  reports,  we  described  the  impact  of  surface  waves  on  low-order  statistical  moments, 
viz.,  mean  winds,  vertical  momentum  fluxes,  and  turbulence  variances  {e.g.,  see  figures  3  and 
4).  In  order  to  link  statistics  with  flow  events  and  structures  we  carried  out  extensive  flow  vi¬ 
sualization  and  constructed  animations  of  select  flow  variables  from  the  database  of  available 
LES  solutions  summarized  in  Table  1.  An  interesting  (typical)  result  is  presented  in  the  three 
panel  plot  of  figure  5  where  we  compare  the  instantaneous  static  pressure  p' / p  and  flow  vectors 
(u  —  uc,  w)  for  simulations  EB,  E,  and  A.  Following  Adrian  et  al.  (2000)  we  choose  the  constant 
velocity  uc  so  as  to  expose  the  rotational  (vortical)  motions  in  the  surface  layer.  In  each  simula¬ 
tion  the  surface  boundary  conditions  imposed  at  the  bottom  boundary  of  the  LES  are  identical, 
viz.,  the  resolved  wave  spectrum  and  the  small  scale  surface  roughness  z0.  All  three  simulations 
are  performed  with  neutral  stratification  (Q*  =  0  at  the  surface)  but  with  varying  geostrophic 
winds  (see  Table  1). 

Inspection  of  the  flow  fields  presented  in  the  x  —  z  vertical  slices  of  figure  5  clearly  shows  the 
impact  of  surface  waves  on  the  turbulent  surface  layer  winds  for  varying  wave  age.  For  the  sta¬ 
tionary  bump  case  EB  (no  moving  waves)  we  observe  that  the  pressure  field  is  strongly  cor¬ 
related  with  the  underlying  wave  shape,  i.e.,  (negative,  positive)  p' / p  is  well  correlated  with 
wave  (crests,  troughs),  respectively,  with  coherent  pressure  signals  extending  upwards  to  20  m  or 
more.  This  is  the  expected  pattern  for  turbulent  flow  over  low  amplitude  hills  in  the  absence  of 
flow  separation  {e.g.,  Gong  et  al.  1996).  The  middle  panel  is  a  typical  flow  pattern  that  occurs 
with  moving  surface  waves  when  the  wave  age  is  near  equilibrium  A  ~  1.5.  In  this  situation  the 
pressure  field  induced  by  the  surface  waves  is  dramatically  weaker  than  in  the  stationary  bump 
case  and  there  is  only  a  faint  correlation  between  p'/p  and  the  wave  shape  h(x,y ).  Essentially, 
most  of  the  long  wavelength  (resolved)  surface  waves  are  traveling  near  the  local  wind  speed 
and  as  a  result  these  long  waves  induce  only  small  perturbations  in  the  surface  winds.  The  ab¬ 
sence  of  a  significant  pressure-wave  slope  correlation  implies  that  the  form  drag  induced  by  long 
wavelength  waves  A  >  (9(10  m)  in  these  simulations  is  small.  In  other  words,  when  winds  and 
waves  are  near  equilibrium  A  ~  1.2  the  bulk  of  the  form  stress  (drag)  is  supported  by  slowly 


lrThe  meeting  papers  Sullivan  et  al.,  (2010a, b)  are  available  at  http://www.mmm.ucar.edu/people/sullivan/. 


propagating  small  scale  waves  {e.g.,  Makin  et  al.  1995).  The  surface-elevated  contours  of  nega¬ 
tive  pressure  surrounded  by  spanwise  rotating  velocity  vectors  suggests  that  the  dominant  co¬ 
herent  structures  in  the  wind-wave  equilibrium  boundary  layer  are  likely  hairpin  vortices  (e.g., 
Adrian  et  al.  2000  and  see  discussion  below).  In  the  bottom  panel  of  figure  5,  a  strong  correla¬ 
tion  between  the  pressure  held  and  the  underlying  wave  shape  re-appears  for  wave  age  A  ~  4.8; 
(negative,  positive)  p' / p  is  now  well  correlated  with  wave  (crests,  troughs),  respectively.  Ani¬ 
mations  clearly  show  these  pressure  patterns  are  connected  to  the  motion  of  the  surface  waves 
as  they  propagate  faster  than  the  overlying  winds.  Wind-aligned  swell  (as  depicted  in  figure  5) 
can  be  dynamically  important  as  it  generates  a  surface  pressure  held  that  reverses  the  sign  of 
the  form  stress.  Similar  to  Sullivan  et  al.  (2008),  on  average  the  pressure  pattern  induced  by  our 
multi-mode  fast  propagating  swell  tilts  slightly  forward  of  the  wave  crests  (i.e.,  into  the  wind) 
resulting  in  positive  form  stress.  This  is  an  indication  of  “wave  pumping”  which  can  lead  to  a 
wave  driven  wind  in  the  neutral  boundary  layer.  Wave  pumping  combined  with  small  amounts 
of  convection  can  further  modify  Monin-Obukhov  surface  layer  scaling  and  the  bulk  boundary 
layer  properties  (see  hgure  4  and  Nilsson  et  al.  2012). 

The  neutral  low-wind  swell  dominated  case,  with  a  broadband  wave  spectrum,  is  a  particularly 
intriguing  regime  with  the  how  dynamics  determined  by  the  competing  effects  of  airflow  turbu¬ 
lence  and  wave  pumping;  the  former  tends  to  destroy  the  organization  induced  by  the  latter.  A 
vast  body  of  literature  (e.g.,  Zhou  et  al.  1999;  Adrian  et  al.  2000;  Christensen  &  Adrian,  2001) 
shows  that  in  shear  generated  hat  wall  boundary  layers  the  dominant  coherent  structures  are 
packets  of  hairpin  vortices.  These  vortices,  which  are  found  using  various  eduction  schemes, 
critically  impact  the  how  dynamics  in  both  the  lower  and  upper  regions  of  hat  wall  boundary 
layers.  We  follow  Adrian  et  al.  (2001)  and  adopt  the  so-called  A „  vortex  detection  method  to 
identify  the  coherent  structures  in  the  low-wind  swell  dominated  MABL.  Mathematically  A c;  is 
the  magnitude  of  the  complex  eigenvalue  of  the  velocity  gradient  tensor  dui/dxj  and  physically 
A d  is  a  measure  of  local  swirling  strength.  In  our  simulations  where  a  coordinate  transforma¬ 
tion  is  used  to  map  physical  space  onto  hat  computational  space  the  eigenvalues  of  the  velocity 
gradient  tensor  are  computed  from  the  expression 

duj  _  (AkdtA 

dxj  d£k  dxj 

where  d£,k/dxj  are  the  metric  components  connecting  physical  and  computational  coordinates. 
The  real  and  complex  eigenvalues  are  found  by  solving  the  cubic  characteristic  equation  implied 
by  (1). 

Figure  6  shows  contours  of  Xci  with  an  overlay  of  how  vectors  for  simulation  A.  The  dark  red 
colors  are  regions  where  the  local  swirling  strength  is  particularly  high.  Also  notice  these  regions 
are  surrounded  by  spanwise  rotating  velocity  vectors  which  adds  confidence  to  our  interpreta¬ 
tion  that  high  values  of  A ci  are  indeed  vortical  cores.  Closer  inspection  of  hgure  6  shows  that  the 
dominate  vortical  motions  rotate  in  a  clockwise  sense,  are  packed  near  the  wavy  surface,  and  are 
primarily  found  in  the  local  wave  troughs,  i.e.,  between  wave  crests.  Animations  further  show 
the  persistent  character  of  this  coherent  structure  and  its  forward  propagation  with  the  waves. 
This  near  surface  vortical  structure  is  distinctly  different  than  the  hairpin  vortex  packets  typical 
of  hat  wall  shear  boundary  layers  and  further  emphasizes  that  MABLs  can  have  unique  dynam¬ 
ics  depending  on  wave  age. 


IMPACT  /APPLICATIONS 


The  computational  tools  developed  and  the  database  of  numerical  solutions  generated  will  aide 
in  the  interpretation  of  the  observations  gathered  during  the  past  HIRES  field  campaign.  In  ad¬ 
dition  idealized  process  studies  performed  with  the  simulations  have  the  potential  to  improve 
our  understanding  of  the  mechanisms  that  lead  to  surface  drag  under  high  wind  conditions. 


TRANSITIONS  &  RELATED  PROJECTS 

The  current  work  is  a  collaborative  effort  between  NCAR,  university  investigators  and  inter¬ 
national  research  laboratories.  Also  the  present  work  has  links  to  the  ONR  DRI  Impact  of  Ty¬ 
phoons  on  the  Ocean  in  the  Pacific  (ITOP). 
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Figure  1:  A  snapshot  of  the  instan¬ 
taneous  wave  field  height  h(x,  y ,  t) 
that  is  imposed  at  the  bottom  of 
the  LES  code,  h  is  built  from  a 
sum  of  linear  plane  waves,  and 
waves  propagate  in  time  from  left 
to  right  according  to  the  linear  dis¬ 
persion  relationship.  The  horizontal 
grid  spacing  matches  the  LES,  i.e., 
Ax  =  Ay  =  2.34  m.  The  color  bar 
in  units  of  meters. 


Figure  2:  An  oblique  3D  cutout  view  of  the  surface  following  computational  grid  used  in  the 
LES  focusing  on  the  surface  layer.  The  extent  of  the  vertical  grid  shown  in  the  x-z  and  y-z 
planes  is  approximately  30  meters.  Note  the  horizontal  gridlines  become  level  surfaces  for  z> 
100  m.  The  surface  mesh  is  overlayed  with  color  contours  of  wave  height  h(x,y)  in  meters. 


Figure  3:  Vertical  profiles  of  wind  speed  (left  panel)  and  turbulence  variances  (right  panels) 
for  different  values  of  wave  age  A  =  Cp/U\q.  Friction  velocity  u *  is  used  for  normalization.  The 
dashed  black  line  is  the  rough  wall  formula  U/u*  =  In (z/z0)/k,  where  k  =  0.4.  Temporal  and 
spatial  averaging  is  used  to  make  the  statistics.  The  spatial  averaging  is  over  horizontal  planes 
in  computational  space,  i.e.,  at  constant  (. 


Figure  4:  Vertical  profiles  of  vertical  velocity  moments  ( w2)/ul  (left  panel)  and  ( w3)/ul  (right 
panel)  for  simulations  over  fast  moving  swell  (red  lines)  and  a  smooth  surface  (black  lines).  The 
vertical  coordinate  (  is  made  dimensionless  by  the  boundary  layer  height  Zi.  Note  the  increased 
velocity  variance  induced  by  moving  surface  waves.  The  smaller  value  of  the  third  order  moment 
is  evidence  that  surface  waves  have  modified  the  vertical  transport  of  turbulent  kinetic  energy  in 
this  sheared-convective  boundary  layer. 


Figure  5:  Vertical  x-z  slices  of  pressure  fluctuation  contours  p'/p  and  velocity  vectors  (u  -  uc, 
w)  in  the  surface  layer  for  turbulent  flow  over  stationary  and  moving  waves:  stationary  bumps 
(upper  panel);  wind-wave  equilibrium  A  ~  1.5  (middle  panel);  and  low-wind  swell  dominated 
with  A  ~  4.8  (lower  panel).  The  constant  velocity  uc  is  chosen  to  expose  the  rotational  motions 
in  the  surface  layer. 


Figure  6:  Contours  of  ACj  and  flow  vectors  (u  -  uc,  w)  in  the  marine  surface  layer  for  a  low 
wind  swell  dominated  situation  with  wave  age  A  ~  4.8.  The  dark  red  contours  of  ACj  are  regions 
of  strong  swirling  motion  (vortices).  Notice  that  the  vortical  motions  are  generally  strongest  in 
the  wave  troughs,  and  the  flow  vectors  show  that  the  vortical  motion  is  clockwise.  The  constant 
velocity  uc  is  chosen  to  expose  the  rotational  motions  in  the  surface  layer. 


